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We extract the Weyl scalars ^I/q and \&4 in the quasi-Kinnersley tetrad by finding initially the 
(gauge-, tetrad-, and background-independent) transverse quasi-Kinnersley frame. This step still 
leaves two undetermined degrees of freedom: the ratio |^oj/|*4|5 and one of the phases (the product 
l^'oj ■ I '54 1 and the sum of the phases are determined by the so-called BB radiation scalar). The 
residual symmetry ("spin/boost") can be removed by gauge fixing of spin coefficients in two steps: 
First, we break the boost symmetry by requiring that p corresponds to a global constant mass 
parameter that equals the ADM mass (or, equivalently in perturbation theory, that poi fi equal their 
values in the no-radiation limits), thus determining the two moduli of the Weyl scalars |5'o|, |^'4|, 
while leaving their phases as yet undetermined. Second, we break the spin symmetry by requiring 
that the ratio tt/t gives the expected polarization state for the gravitational waves, thus determining 
the phases. Our method of gauge fixing — specifically its second step — is appropriate for cases for 
which the Weyl curvature is purely electric. Applying this method to Misner and Brill-Lindquist 
data, we explicitly find the Weyl scalars vJ/q and \l/4 perturbatively in the quasi-Kinnersley tetrad. 

PACS numbers: 04.25.Dm, 04.30.Nk, 04.70.Bw 

I. INTRODUCTION AND SUMMARY 

Theoretical understanding of sources for gravitational waves and their emitted waveforms is an important step 
for their detection and for gravitational-wave astronomy. Such theoretical studies are typically done as numerical 
simulations. A numerical relativity simulation comprises of three major steps: First, one needs to construct the initial 
data for the physical system of interest; Second, these initial data are evolved in time jlj, and third, the evolved 
data are to be interpreted, and the outgoing waveforms are to be extracted. This paper addresses the last step, 
specifically the problem of wave extraction. It is, however, also closely related to the first step, namely the problem 
of initial-data determination. Indeed, although our method is applied in this paper to initial data sets (as they are 
available analytically — for the specific problem of interest here), it is equally applicable also for evolution data. This 
paper is a continuation of its prequels [3, |3, |3j S @| , and is based on the method developed therein. 

The long-term goal of numerical relativity is to evolve a binary system over many orbits, starting with slow, adiabatic 
inspiral, going through a plunge, and culminating in a ringdown of the produced black hole. Despite dramatic recent 
advances [7], it is unlikely that current numerical relativity codes will be able to to evolve hundreds of orbits in the 
very near future. The alternative is to start the simulation with the binary being at small orbital separation, so 
that only few orbits are left before the plunge. This alternative presents the construction of initial-data sets with a 
fundamental problem. Specifically, it is unclear how to encode in the initial data the gravitational waves that would 
emerge from inspiral from large orbital separations in the past of the initial time for the simulation. In fact, different 
constructions of initial data for the same physical system lead to distinct physical solutions, with different future 
evolutions and different gravitational waveforms [l], Q . The difference between such initial-data sets is basically the 
amount of spurious gravitational waves that are present at the initial time. Lack of control over how much "junk" 
gravitational waves are present at the initial time implies that even if the numerical evolution were done sufficiently 
accurately, it would typically be unknown which physical system that solution represents, if the "junk" data are strong 
enough. 

The last step of the numerical relativity simulation requires the extraction of the outgoing waveform. The latter — to 
coincide with the waves that detectors such as LIGO or LISA would measure — needs to be extracted in the Kinnersley 
tetrad. If the Kinnersley frame is not identified, the extracted Weyl scalar ^'4 would generally be a linear combination 
of all five Weyl scalars of the Kinnersley frame, and therefore include not just the information about the outgoing 
transverse radiative degrees of freedom, but also information about incoming modes, Coulombic data, and pure gauge 
(longitudinal) information. Indeed, it was shown recently that the extracted waveform changes markedly when the 
tetrad basis vectors are transformed @| (see also Q ) , so that if the Weyl scalars are extracted in a tetrad that is not 
sufficiently close to the Kinnersley one, the predicted waveform would deviate significantly from the waves detected 
experimentally. Identifying — from numerical relativity data — the physically relevant tetrad in which the Weyl scalars 
are to be extracted is therefore an important problem, and this papers — like its prequels — addresses it. 

An important element in the proposed method for wave extraction is the quasi-Kinnersely (qK) frame. A frame 
here means an equivalence class of tetrads ( "particular set of four basis vectors" ) that are related to each other by type 
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III ("spin/boost") rotations (and exchange) transformations ^^l. The quasi-Kinnerley frame is a frame that evolves 
continuously into the asymptotic Kinnersley frame, as the simulated (generically Petrov type-I) spacetime settles 
down to quiescence (generically Petrov type-D spacetime). A particular subset of all quasi-Kinnersley frames is the 
transverse quasi-Kinnersley frame, defined by the requirement that it is a quasi-Kinnersley frame for which the Weyl 
scalars and ^"3 vanish. Transverse frames come in threefold in type-I spacetimes, the transverse quasi-Kinnersley 
frame being one of these three. The transverse quasi-Kinnersley frame is unique for any generic (type-I) spacetime, 
and hence lies its importance. In Ref. [3| we showed how the transverse quasi-Kinnersley frame can be found from 
spatial data. The analysis of Ref. Q still leaves an unbroken residual symmetry, namely the (continuous) spin/boost 
rotations (and the (discrete) exchange operation). In this paper we address the problem of breaking this residual 
symmetry. We do not develop a general method in this paper, that would parallel the approach of Ref. [sj. Instead, 
we take the approach of addressing a particular (simple, but not too simple) spacetime, and attempting to ad hoc 
break the residual symmetry. We believe that this paper sheds light on the general symmetry breaking problem, as 
we discuss below. 

Specifically, in this paper we consider two initial data sets for the time-symmetric (momentarily stationary) head- 
on collision of two non-rotating black holes, namely Misner [l3| and Brill-Lindquist TT] data. These initial-data 
sets differ in the topological properties of the geometry extended through the black holes' throats. For an external 
observer, the two initial-data sets differ in the amount of distortion that each black hole suffers because of the presence 
of the other black hole. Time symmetric black hole binaries were popular in the early days of numerical relativity, 
and their gravitational-wave contents have been compared using other methods 12, 13]. In this paper, we revisit 
this question using a new method, that we believe to be robust: some questions get answers that are background-, 
coordinate-, and tetrad-independent, in addition to being fully calculable from pure spatial data. Other questions 
requires gauge fixing in order to be answered, but this gauge fixing can be done in a robust and generic way for at 
least some problems. 

The benefit of time-symmetric black hole binary data is that they allow for analytic (perturbative) solutions, at least 
for small black hole separations. Specifically, the magnetic part of the Weyl curvature vanishes for time-symmetric 
data. Not only does this property simplify the derivation of the solution considerably, it also facilitates the gauge 
fixing that is required in order to find the quasi-Kinnersley tetrad. Unlike earlier works, we can answer the question 
of which initial data set has more "junk waves" without actually evolving the data: our method allows us to find the 
energy flux in gravitational waves using pure spatial data. 

For either initial data set, we obtain the solution for the Weyl scalars ^'o and ^'4 using a combination of gauge 
invariant and gauge fixing techniques. Specifically, we first find the transverse quasi-Kinnersley frame, which is a 
family of (null) transverse tetrads connected by a spin/boost rotation. The transverse quasi-Kinnersley frame can be 
found in a completely gauge-, tetrad-, and background-independent method from pure spatial data as are commonly 
used in numerical relativity. The residual spin/boost symmetry implies that there are four degrees of freedom to be 
determined: the two moduli and the two phases of the Weyl scalars and ^'4. However, already at this gauge- 
invariant level two of these four degrees of freedom are determined: we determine the product of the two moduli and 
the sum of the two phases by finding the so-called Beetle-Burko (BB) radiation scalar ^ [2], which is invariant under 
spin/boost rotations. Recall, that in a transverse frame ^ = ^'4, so that the phase of the BB scalar ^ equals the 
sum of the phases of the two Weyl scalars, and its modulus equals the product of the moduli of and 5*4, when the 
latter are in a transverse (quasi-Kinnersley) frame. 

We break the remaining spin/boost symmetry in two steps: First, we fix the boost degree of freedom by requiring 
that the expansion (or spin coefRcient p) corresponds to a global mass parameter that equals the system's ADM 
mass [Til . This method is equivalent for the specific problem of interest to the requirement that p (or, equivalently 
fi) coincides with the no radiation limit (in which the two black holes degenerate into one). While this gauge fixing 
is problem adapted, we believe it — or an adaptation thereof — is general enough for most applications in numerical 
relativity simulations. Specifically, we require that the spin coefficient p (or, equivalently, the expansion 9) coincides 
with that of a Schwarzschild black hole when the proper distance along the spacelike geodesic connecting the throats 
of the two black holes vanishes. This method is applicable also to spacetimes that include perturbed rotating black 
holes. This gauge fixing determines the two moduli |^o| and |5'4| completely. 

The remaining degree of freedom can also be determined by gauge fixing, at least for the specific problem of interest 
here, although not without additional information. Specifically, we argue that the spin rotation can be fixed by 
choosing the ratio of the spin coefficients tt (which below we denote by w) and r. In principle, each of these spin 
coefficients has a non-zero spin weight and vanishing boost weight, and could be used to fix the remaining degree 
of freedom. However, both have zeros, where spin rotations leave them vanishing. Consequently, gauge fixing is 
not unambiguous. However, their ratio exists everywhere (for the initial data sets of interest in this paper), and 
has vanishing boost weight and non-zero spin weight, so that its use is vary natural, and completely determines the 
Weyl scalars 4*0 and '^4. Specifically, we extend results from point particle motion and expect that only the /i+ 
polarization state is excited significantly. Indeed, most studies of point particle motion near black holes show that 
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the hx polarization state is extremely hard to excite [15[. As the motion of the two black holes in the problem of 
interest is radial, indeed no excitation of the polarization state is expected. We therefore choose the gauge fixing 
so that the generated waves include the /i+ polarization state only. 

Our proposed gauge-fixing requires knowledge of the polarization state of the emitted gravitational waves. While 
the total emitted energy in gravitational waves is independent of this second step of gauge fixing (and therefore the 
flux of energy can be found without first breaking it), its distribution between the h+ and hx polarization states is 
not. We argue that a physical assumption on the polarization of the waves needs to be made. Our method of gauge 
fixing, specifically its second step, is not specific to the head-on collision case. We expect it to be appropriate for 
cases for which the Weyl curvature is purely electric. When the magnetic part does not vanish, tt and r may have 
non-coinciding zeros, thus hindering the use of their ratio. The breaking of the residual spin symmetry when the Weyl 
curvature has a non-vanishing magnetic part awaits further consideration. 

The organization of this Paper is as follows: In Section[TT]we review the construction of the quasi-Kinnersley frame. 
We then construct this frame — and break first the boost and then the spin symmetry to find the quasi-Kinnersley 
tetrad — for Misner (Section llTl|) and Brill-Lindquist (Section [rv]) initial data sets. Finally, we compare our results 
for Misner and Brill-Lindquist initial data in Section [Vl 



II. CONSTRUCTION OF THE QUASI-KINNERSLEY FRAME 

In this Section we outline the calculation of the quasi-Kinnersley frame from spatial data. See Ref. [s*] for more 
detail. In a 3 -|- 1 decomposition, the gravitational field is expressed in terms of the spatial metric 7^ and the extrinsic 
curvature K^j . From these, the electric and magnetic components Eij and Bij of the Weyl tensor can be computed 
from 



£^J = R^J+KK^J-K^kKJ''-4^TS^J (1) 
E^, = - ^l^Jl'''Skl (2) 

where Sij = ^^"'^j'^Tab is the spatial projection of the stress energy tensor Tab, and 

= ^t^'VkKi, (3) 
= (4) 

Here, e^fe := T'^emijk is the intrinsic spatial volume element, where r™ is the unit future-directed normal to the 
hypersurface, and emijk is the spacetime volume element. 

By construction, Eij and Bij are both symmetric and traceless. We then define the complex tensor 

(5) 

from which the scalar curvature invariants / and J can be computed as 

1 

-( 

2 
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I = 7;C',C', (6) 



J = —C^^C\C\,. (7) 



6 ~ ' ' 

In terms of the dimensionless Baker-Campanelli speciality index ( "the S invariant" ) [if 



S = 27^ (8) 



the Coulomb scalar can be written as 



o,±_ 3JW^(5)V3 + w^(5)-V3 



u,± ^ __ xv-y ' --xv-y /g^ 



where W^{S) ^ \fS - VS - 1, and the BB scalar i 
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^o.± = i / [2 - W^iSy/"^ - M^5(5)~i/3] , (10) 
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where W^{S) = 25' — 1 + 2^/3(3 — 1). In general, both the Coulomb and the BB scalars admit three distinct complex 
roots, so that ^ and x are multiply defined. We denote these three roots with the superscripts "0, ±" . This ambiguity 
is related to the fact that generic (i.e., algebraically general) spacetimes admit three distinct transverse frames, in 
which the Weyl scalars take different values. The value associated with the principal branch will be denoted with 
a superscript 0. The other two branches with superscripts ±, respectively. The quasi-Kinnersley 4] value will have 
no superscript. For the cases studied here, there is a clear limit in which spacetime is algebraically special (more 
specifically, is Petrov type-D). One can then find the eigenvector CTq corresponding to the quasi-Kinnersley frame 
and the BB scalar in the quasi-Kinnersley frame using the methods of Refs. [1, @. Specifically, the quasi-Kinnersley 
frame can be identified by choosing the eigenvector ag of corresponding to the eigenvalue of greatest modulus [1] . 
Notably, the BB scalar is not needed explicitly in order to find the quasi-Kinnesley frame. The only required quantity 
is the eigenvector ctq . We emphasize that while the BB scalar is a promising tool to (partially) describe the radiative 
degrees of freedom, the wave extraction method is not based on finding it. Rather, it is an auxiliary quantity, which 
need not be found for determining the quasi-Kinnersley frame. 

Once the eigenvector has been found, one might like to find the elements of some tetrad in the quasi-Kinnersley 
frame. This, too, can be done using only spatial, rather than space-time, data. The real null vectors of an arbitrary 
tetrad project into a space-like hypersurface to define a pair of real spatial vectors. While the normalizations of these 
spatial projections naturally vary if one performs a spin/boost on the tetrad, they define a pair of invariant rays in 
the tangent space at a point. These rays correspond to a pair of real unit vectors A" and z)" which are generally 
completely independent of one another. In particular, they point in exactly opposite directions only when the normal 
t'^ to the spatial slice lies in the space-time tangent 2-plane spanned by and n'^. This degenerate case, which 
happens when the magnetic part of the Weyl tensor vanishes, is treated separately in Ref. 's'l, and is applied below 
for the moment of time symmetry of Misner and Brill-Lindquist initial data sets. The fiducial time foliation used in 
numerical relativity is overwhelmingly likely be "boosted" relative to the 2-plane defined by a given null tetrad, and 
only parallel A" and i)" are actually disallowed. 

Since A" and j)" determine the directions of the space-time vectors i"" and n'^, they also suffice to determine the 
frame associated to (T°, and thus a" itself. The explicit relation is 

a" = (1 - A • t>)-i (i>" - A"^ + ie"'"' Xt v^). (11) 

Note that when A"^ and v'^ are interchanged, ct" simply changes sign^s it should since the bivector Eaf, changes sign 
under an exchange operation in the corresponding frame. (See Ref. [3| for definitions and detail.) 
To invert Eq. (fTTj) and solve for A" and v"^ given (T°, we separate ct" into real and imaginary parts: 

a" ^ x'' + iy" . (12) 

In general, neither nor is unit, but the normalization condition for ct" demands they be orthogonal with norms 
satisfying — ~ 1. Taking the electric part of Eq. (10) of Ref. [3] then yields 

,a&c a abc i a 

= NP ' ^ ■ ^ ^ 

Finally, the elements of the most general null tetrad (i?°, n", to", to") in the frame associated to ct" take the form 

r = j'^l'! (f" + A") 



- _l£L 



(f"-hz>") (14) 



where the complex unit projection of to" is 

A" = (1 + A • i))-^ (A" + i>" + ie"''" h Vc)- (15) 

These vectors are the result of a spin/boost with parameter c — |c|e*'', acting on a preferred tetrad fixed by setting 
^" fa = n° fa < and m" fa = m° fa < 0. Note that when A" and j>" are indeed anti-parallel, the expression for to" in 
Eq. (fT4|) is ill-defined. However, the limit as i>" —A" does exist, though it does depend on how the limit is taken. 
This is natural since the fixed tetrad is determined only up to spin transformations in this case. This degenerate case 
is treated in detail in Ref. @, and is used below for the head-on collision case of two nonrotating black holes as an 
illustration. Notice that Eqs. (jl4p are determined up to a continuous complex function c, that is, the basis vectors are 
a two (real) parameter family ("equivalence class") of tetrads. This family of tetrads is the quasi-Kinnersley frame, 
and the desired quasi-Kinnersley tetrad is a single member of this family of tetrads. 
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III. MISNER DATA 

Misner data can be written in spherical isotropic coordinates as 

da^ = ^^{R, e- Ho) [ dR^ + R^ ( dO^ + sin^ 9 d(j)^)] (16) 
with unit lapse and vanishing shift, where the conformal factor is given by 

i+i 

R 



$Misner(i?,^;/X0) = 1 + 2 ^ ( ^ ) Pt{cOs9), (17) 



e=0,2A,--- 

where P£(cos^) are the Legendre functions. The parameter fiQ is related to the separation of the two holes (in a 
manner to be made clear below) and M is the ADM mass of the system. The coefficients Ki{fio) are given by 

jj^t ^-^^^^^ and ^.-.^t-^. (18) 
(42^1 '^-^ smhn/in smhniin 

^ ' n—\ n— 1 

As pointed out by Anninos et al. (T^ , Misner data can represent perturbations over a single Schwarzschild black hole 
for Ho ^ 1. Repeating the approximations in Ref. ^l5|, for £ > 1, we assume that (cothn/io)^/ sinhn/iQ ~ {nfioY 
(which strictly speaking is correct for n/zo <C 1) for small /io, so that we only keep terms in summations up to ~ A*o ^■ 
We therefore have 

^ (cothr^ _ C(i+1) ^ llHA^ol , . .^Cii+l) 

> ^-r ~ — jTT- ' ^1 ~ -and Keipo) ^ — . 19) 

^ smhn^o |41n^or+ 

Here, C is the Riemann zeta function. Note that Si = M/2. 

The proper distance L along the spacelike geodesic connecting the throats of the two black holes is 



L^2 



1 + y — 

■"^^ sir 



^ sinh(n/^o) 

n— 1 

which in the limit of small separation {/iq <C 1) is given by 



(20) 



L«— -M, (21) 

4| In^ol 



so that the initial data, to hexadecapole order, can be written in this approximation as 

M C(3)P2(cos6') flV C(5)P4(cos6i) f ^ 



" ^+2i^+ 2.^ [r) + 2.10 [r] ' (22) 

A. The quasi-Kinnersley frame 

We first find the elements of the tensor C^j to 0{L^). They are 

^ _ 7 MR^ 2H(3) p. ... 2AR^ + 2MR + M' 

- WTW'^^ ^ — {2R + My — ^ +^(^ ) (^^) 

_ .6 MR^ , 25C(3) 2P2(cosg)(14i?^ + 3Mi? + Af^)-(2fl + M)^ ^3 , 
^0 - ^ W+Mf ~ {2R + My ^ +^(^ ^ ^^^> 

^0 _ .6 , 2^C(3) 2P2[cose){lQR^ - MR) + {2R + Mf , 



_ 1 r^R 
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The eigenvalues of C*^ are 

Ai = + 0{L^) 
A2 = + 0(l6) 

A3 = 

so that the eigenvalue that corresponds to the quasi-Kinnersley frame is Ai, and the corresponding eigenvector is 
v1 =a{S° + Xm S^), where 



X,r = sine C0S9L' + 0{L^) 



and 



27r6 Mi?4 



8C(3) 



(2i? + M)2 7r6 (2i? + Af )3 



P2(cos6') L3 + 0(i^). 



(27) 



(28) 



(Notice, that Xm affects a only at 0{L^).) 
Following the method outlined in Refs. [1, [1], we find the BB scalar to equal 



M 



2«-32 Cn3) 



sin"* 61 



5.29.3C(3)C(5) 



7ri2 {M + 2i?)io 
and the Coulomb scalar to equal 

2^ M 2^ C(3) (24 i?^ + 2 Af i? + A-P 



Trie i?2 (M + 2R) 
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[7 P2 (cos 61) + 2] sin'' 6*2.*^ + 0{L^) 



(29) 



XM 



(M + 2Rf 



7r6 (Af + 2i?)7 



3C(5)(10i?2 + 5A/i? + A/2) 
ttIO i?2 (Af + 2i?)7 



(30) 

Because this is a time symmetric configuration, the magnetic part of the Weyl tensor vanishes, from which it follows 
that the eigenvector of Weyl is real, as is evident from Eq. (I27p . This implies that the two real projections of the two 
real null vectors of the quasi-Kinnersley frame onto the hypersurface are anti-parallel. As v" is already normalized, 
we take A" = — i>" = v^, and then the two null spacetime vector are 



V2 



1 



(s^+xms^: 



(31) 



V2 



6^ - 



and the complex null vector is 



m 



V2$2^ 



smy - yrrw^ 



(J^ + ^M^e") 



Here, c := |c| exp(ii?) is a complex valued function on spacetime with modulus |c| and argument 'd. 
The Weyl scalars in the basis (I5T]) - ([55|) are given by 



2 2,;^ 3 • 24C(3) sin2 



c -e 



|cr2e2^^ 



(2i? + M)' 



-2^2.. C(5) ,ij,2^[7p^(^^^^) 



L5 



R^(2R + Mf 



2 r(3) . 2a2P2(cos6')(64i?2 + 49Mi?-|-A//2)-|-(2i?-hAf)(32i? + Af) 



sin^ 61 



MR^{2R + Mf 



L^ + 0{L'^ 



2^ Af i?3 2^ C(3) (24 p2 + 2 MR + AP 



(Af + 2Ry 



TT^ (Af + 2Ry 
,2 /) 



P2(cos6l)L-^ + 0(L^ 



l^|2g-2,^ 3 ■ 2 C(3) ^^^^^ ^3 ^ |^|2g-2.^ sin^ ^ [7P2(cos0) + 2] 



(2P-H Af)5 



rlO 



£5 

P2(2P + Af)5 



+ ere 



2 -2,^ 2'^C (3) 2 . 2P2(cos6i)(64p2 + 49Mi? + M2) _^ Af)(32P + Af) 



MR3{2R + MY 



(32) 



(33) 



(34) 
(35) 



P*^ + O(P^) . (36) 
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Because this is an approximate solution, the longitudinal Weyl scalars ^'i and ^3 are not exactly zero, but at O(L^). 
Specifically, as we expressed Xm to 0{L^), the coefRcients of and \E'3 at 0{L^) vanishes, and the next term, at 
0{L^) remains. We can nullify these (and successive) terms by keeping terms to sufficiently high order on $ and in 
Xm- 



B. The qK tetrad: Breaking the boost symmetry 



We break the spin/boost symmetry in two steps: first, we break the boost symmetry while leaving the spin symmetry 
unbroken. At the next step we address the breaking of the spin symmetry. To break the boost symmetry we are 
looking for a spin coefficient with non-zero boost weight but vanishing spin weight. The two candidates are the spin 
coefficients p and fi. The spin coefficient 



= -2V2 



1 R{2R 



M) 



{2R + Mf 



0{L^). 



(37) 



Under Type III rotations p transforms like i.e., p \c\~ p. We can therefore scale the null basis vectors ^" and 
n" by demanding that p equals, in the limit L = when no radiation exists, the Schwarzschild black hole value ps- 
As p is real, the twist vanishes, and p equals the expansion 6. Recall that for L = 0, we have a single non-rotating 
black hole with ADM mass M, for which the spin coefficient ps = — i?~^(l + M/2i?)~^. Requiring that for L = we 
get p = Ps, we find that 



1 2R-M 



sj2 2R + M' 



(38) 



Alternatively, we can recall that for L = (also for a Kerr black hole), '^2/ p^ = M. Under a type III rotation the 
combination '^2/p^ has boost weight 3 and zero spin weight. We can therefore find a boost rotation function |c| so that 
this combination equals a global constant. To fix the remaining multiplicative constant factor in |c| we require that 
the global constant equals M. (In practice we apply this method here in the L ^ limit.) Putting these expressions 
together, we solve for |c|, and recover the previous result. One can also use the spin coefficient p := ffi^m^ V/jriQ,, and 
demand that it would equal the no-radiation limit as L ^ 0, or ps = —2R{2R — M)^/{2R + M)"^. This calculation 
recovers the value of |c|. 

We can now find 1^*0 1 and |^4| in the quasi-Kinnersley tetrad 



3 • 2^C{3) 



2^2(3) 



rl2 



sin^6l 



{2R - Mf{2R + Mf 
2P2{cos 9) {64R^ +49 MR 



+ 



5 • 25C(5) 



1 



M' 
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R^{2R- Mf{2R + Mf 



sin2 6'[7P2(cos6') -F2] 



MR^{2R + M)4(2i? - M) 



{2R + M){32R + M) ^ 



(39) 



l*4l = 



3-23^(3) ^.^2 g {2R-MfL^ ^ f>-2^(^{b) {2R- M) 



{2R + My 

2'C^(3) „;n2 2P2(cosg)(64i?2 + 49MR + M^) + {2R + M){32R + M) M\^T^^n(T^\ 
- ^^^1° ^ Mi?3(2i? + M)8 {2R-M) L + 0{L ) 



(40) 



At great distances R:^ M, |*4| ~ L^/R^ and |*o| ~ L^/R^. 

The energy flux in gravitational waves depends on the (time domain) Weyl scalars \I/o and \l/4 in a non-local way, 
as is manifest from the integral over preceding times in the expression for the energy flux in gravitational waves (at 
infinity) 



dE , , 



■ lim 

r^oo 



dt'ip{t',r,9,(l)) 



(41) 



where ip is the Teukolsky function, r is the geometrical radial coordinate, and du) is the integration measure over 
the 2-sphere. The contribution to the energy fiux (at finite distances) in gravitational waves at the moment of time 
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symmetry ( "t = 0" ) from that moment in time ( H' = 0" ) is given by 



^0 := j R'\^o?d^l 

3-2"C^(3) 3-2"C^(3) MR^ + 22MR + 



Stt" R^{2R + Mf{2R- my S-TttI^ MR^{2R + Mf{2R - M) 



5-2"C^(5) 2i^C^(3)C(5) 160i?+173M „ 

7ri9 R^{2R + Mf{2R- MY 5 • 77r2i MR^{2R + Mf{2R - MY ^ > ^ > 

Ti := I R^\-^A^dn 



3-2^(3) {2R~MY 6 3-2^C^(3) (64i?2 + 22Af + M2)(2fi - M)^ „ 

Stt" i?2(2i? + M)io 5-77ri7 Mi?5(2i? + M)" 

5-2V(5) (2i?-Af)4 p 2V(3)C(5) (160i?+173Af)(2i?-Af)4 

7ri9 i?6(2i? + M)io 5-77r2i MR%2R + My^ + I J' I i 

where JF^ is the flux in 

C. The qK tetrad: Breaking the spin symmetry 

To break the spin symmetry we look for spin coefficients that have non-zero spin weight and vanishing boost weight. 
The two candidates are the spin coefficients zu^ and r. The spin coefficient 

vj := m''tVi,na 

and 

We cannot just consider a spin rotation (i.e., a type III rotation with |c| = I) directly on tu and r, because -cc and 
T have zeros on the equatorial plane and at the poles. At the zeros (and numerically also in neighborhoods thereof) 
the rotation parameter is arbitrary. However, because both w and r have zeros at the same locations, the limit of 
their ratio exists everywhere. We therefore consider the ratio 

- = 6-2^" . (46) 
r 

Notably, Eq. (HB)) is satisfied nonperturbatively. Under type III rotation by angle zu/t ^ e~'^'^^zu/T (i.e., the ratio 
w/t has spin weight —2), so that by choosing a spin rotation function ip we can determine the ratio w/t uniquely. 
Choosing the function ip we in fact determine the local polarization state of the gravitational waves in the wave zone, 
where ^'4 and ^Pq have a clear physical interpretation in term of outgoing and incoming waves, respectively. Notice, 
that once the boost degree of freedom is fixed, the total energy flux in gravitational waves is determined. The spin 
degree of freedom will nevertheless move part of this flux between the two polarization states. This property suggests 
to us that the spin degree of freedom can only be determined based on the understanding of the polarization of the 
waves in the wave zone. As already discussed in the Introduction, studies on radial motion in the point particle limit 
have shown that the /ix polarization state is extremely hard to excite. Following Berti et al [3| we may presumably 
lean on the point particle results, and assume that only the polarization state is generated. In our case, indeed 
this conclusion is guaranteed by the diagonality of the metric p^. 



^ Conventionally, m°'£^ Vj, ria is denoted by n. Here, to avoid confusion with the geometrical n, we change the usual notation, and denote 
it with ro. 
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Our ability to break the spin symmetry depends on the observation that the ratio w/t is free of zeros, which may 
be attributed to the vanishing of B"^^.^ We assert that this method for breaking the spin symmetry by gauge fixing 
w/t according to the polarization state is applicable to purely electric cases for which = 0. The question of how 
to break the spin symmetry for the magnetic case is not answered by the above discussion, and we hope to address 
this question elsewhere. 

In practice, both and 5*4 would be real for this case. Choosing = —d we can set tu/r = 1. Under this rotation, 
"Ho e'^^'f^o and *4 -> e^^i^^^ ^ so that the new Weyl scalars equal their moduli, i.e., 

^ 3 ■ 3) 2 g I ^ • ^y 5) 1 -sm^e[7P2icos9) + 2]L^ 

2"C^(3) 2n '^P2{cosemR^ + 49MR + M^) + {2R + M)(32R + M) , 
- ^ Mi?3(2i? + M)4(2i?-M)2 ^ +^(^ ) ^4^) 



*4 = 



3 • 23C(3) g, (2^ - ' 5-23C(5) {2R~ M) 



{2R + My 



rlO 



22(2(3) . 2,2P2(cose')(64i?2+49Afi?+M2 
77; — sm 



-477 -T^ sin^ e [7P2(cos 6*) + 21 L 

R^{2R + My ^ ^ 

{2R + M){32R + M) 



MR^{2R + Mf 



{2R - Mf + O(L^) . (48) 



We emphasize that this breaking of the spin symmetry relies on physical intuition based on extension of the point 
particle results to the equal mass case. Any other choice for Lp would do just as fine in the absence of knowledge of the 
local polarization state of the gravitational waves. Notice, in contrast, that breaking the boost symmetry is a unique 
gauge fixing by virtue of the radiation free limit, and that getting the quasi-Kinnnersley frame is gauge invariant. 
Finally, we can write explicitly the basis vectors for the quasi-Kinnersley tetrad: 



2R + M 
2R-M 



S? 



1 



(49) 



2R-M 
2{2R + M) 



{51 + Xm5'S) 



(50) 



m 



1 



V2$2^ 



(51) 



V2$2i? 



■51- 



^S-g-R'XMS'k) 



(52) 



IV. BRILL-LINDQUIST DATA 

Next, we write Brill-Lindquist initial data. The metric has the form (|16p . but now the conformal factor is given by 



$BL(i?,^;^o) = l + y 



' i?2 sin2 9+{Rcose~ zo)2 J R^ sin^ + {R cos0 + zq 



(53) 



where the two black holes are at positions ±zo along the z axis. The parameter m is related to the mass of the black 
holes. 



^ Where B'- ^ the real unit vectors A° and are not co-linear, so that the limit ro/r may not exist everywhere. 
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A. The qucisi-Kinnersley frame 

In this Section we find the quasi-Kinnersley frame for B-L data. Our calculation parallels its Misner data counter- 
part, so we describe it briefly. 

= -^WTW^"^ (R + mr P^i--^(^)4 + Oi4) (54) 

= ' (R + Lr ^ ^^^^ + ^"^^ + ^;;y^jr'^'^''^"^' ^'0 + 0(4) (55) 
, mR^ R{5R - m) P2(cosg) + (i? + mf 

^ - ^ii + m)6 + ^™ {R + my ^° + ^^''> 

m + 4i? 

= "^"^ E(i? + m)6 '^"^ ^ ^ ^0 + ^^^0 ^ ^^^^ 
The eigenvalues of C are 

M=C\ + 0{zt) 
X2 = C% + 0{zt) 

so that the eigenvalue that corresponds to the quasi-Kinnersley frame is Ai, and the corresponding eigenvector is 
=a(5« +XBi5^), where 

Xbl = sin e cos e zl + 0{z^o) (58) 

and 

a = ^-(1 + R^Xl,)-^'^ = - 2 P.(cos.) zl + O(.o^) . (59) 

(Notice, that Xbl afl'ects a only at 0{zq).) 
We repeat the calculation for Brill-Lindquist initial data, and find that the BE scalar equals 

^BL = rj?!^\^o sin' zl + ^3.^""' [(3 - 14 m i? - 2 m2) P2 (cos Q) - {%R^ mR^m^)\ sin^ ^ 4+0(^,^) 

(60) 

and the Coulomb scalar is given by 

2mi?3 2m(6R'^ + Rm + m'^) ^ ^ 9 4x 

= -JrT^ -iRT^r ^^(cos 0) zl + O(.o') . (61) 

The Weyl scalars in the quasi-Kinnersley frame are 
*o = Icp^e^^'' IS^^^ zl + \c\-^e'''^ .""J^"'^,. [(SP^ - 14mP - 2m^) P^icosO) + {6R^ + 7mR + m^)] 4 + 0(^) 



= -'WT^ -'"^ (R + my ^^(^os.).g + Q(4) (63) 
*4 = |c|^e-2^'' ^^^^ zl + \c\'e-'^ \"''l'\, m' - 14mP - 2m') P^icosO) + {6R' + 7mR + m')] 4 + 0(S) 
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B. The quasi-Kineersley tetrad: Breaking the boost symmetry 



For zo = 0, ps = —R/{R + mf'. Requiring that the spin coefBcient p equals its counterpart in the no-radiation 
Umit, we find that 



1 R-m 



V2 R + m' 



(65) 



so that in the quasi-Kinnersley tetrad 



l*nl = 



6m sin 9 



2m sin 9 



{R + mY{R-mf ° R^{R + mY{R-mY 
+ 0{zl)\ 



[(3i?2 - 14mi? - 2m^) P2(cos 61) + {QR^ + 7mR + m^)] 

(66) 



1*4 



3m(R — m)^sin'^6 ^ m{R~ siv? 9 „o „ os „ , , „o „ on-, a 

^ ' ~2 , V ) [{m^ -lAmR-2m^) P2{cos 9) + {&R^ + 7mR + m^)]z^ 



2{R + my 



^° ^ 2i?3(i? + m)8 



(67) 



Again, we could also determine the boost degree of freedom using the combination '^2/ P^- 

Finally, the contribution to the energy fiux in gravitational waves at the moment of time symmetry from that 
moment in time is given by 



^4 



:= j R^\^, 



3 ■ 



■ Zn :: — - / ^ no / ^ — Zq -|- U\Zq) 



f>R^{R + mY{R-mY ° b-lR^{R + mf{R-mY 



:= j i?'|*4 



3-2Vm2(i?-m)^ 4 3 • 2^Trm^{R - m)^{16R^ + llmi? + m^) g s\ 



5i?2(i? + m) 



5- 7i?5(ii + m)ii 



(68) 



(69) 



and 



C. The quasi-Kineersley tetrad: Breaking the spin symmetry 

The spin coefficients 

n7 = -V2zsin6>cos6' „,/„ , r z^ + 0{z^) (70) 

it^(i£ + m) 

r = -^/2ism9cos9^^^^^zl + 0{z^) (71) 

so that again we can find the ratio (whose limit exists also on the equatorial plane and at the poles) to be 

- = e-^'^ . (72) 

T 

Doing a spin rotation with, as above, parameter ip = —-d we can set the ratio vc/t to equal unity as the gauge fixing 
that gives only h+ gravitational waves, and does not excite the hx polarization state. In this gauge, the Weyl scalars 
are given by 



6m sin^ 9 2 
° " (i? + m)3(i?-m)2 ^0 



2?77 SlTl^ ^ 

+ -537^ ^4r5 ^ [(3^' " l^mi? - 2m^) P2{cos9) + {6R^ + 7mR + m^)] 4 + 0(^0') (73) 

R'^{R + mp{R — m)'' 



_ 3m{R - to)2 sin^ 9 2 
2{R + my ^° 

+ "'2R%R + my^ t^^^' ~ ^^""-^ ~ ^""'^ + + ^""^ + ""'^^ + ■ ^^^^ 
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V. COMPARISON OF MISNER AND BRILL-LINDQUIST DATA 



We cannot yet compare our results for the Weyl scalars and ^'4, or even our results for the BB scalar for the 
two cases: the expansion parameter zq for the case of Brill-Lindquist initial data is different from the proper distance 
L we used in the Misner case; we need to relate Zq to the proper distance between the two holes. One way to do 
this is to find the proper distance between the apparent horizons of the two holes. Alternatively, we require that the 
Coulomb curvature is the same for both initial data sets (to quadrupole order). Comparing xm and xbl to second 
order, we find that the relations between the parameters of the two sets are: 



M = 2m and ^ = m z?. . (75) 
2 7r6 " ^ ' 



Using these two relations, we re-write 



and 



e„ - 9"^' sin^ ez^ + 5 2^/3 3 ^^^^ (TP^ (cos g) + 2) , 



ri2(M + 2i?)io 

^ ^ ^ [2{6R^ -UMR- M'^)P2{cos9) - {2AR^ + UMR + M^)] sin^ 6 + 0{L^^){77) 



We see, that the difference between the BB scalars for the two cases appears at 0(L°), and arises from the product 
of the 0{L^) of one Weyl scalar with the term at 0{L^) of the other. One may therefore expect the energy flux in 
gravitational waves of the two cases to also differ at 0{L^). However, we find that they differ at 0{L^'^). The reason 
is that the contribution of the 0{L^) term to the energy flux in gravitational waves vanishes after integration over the 
sphere. The next term, at O(L^), arising from the product of the 0{L^) and 0{L^) terms, is equal for both initial 
data sets. At 0(L^°), which has contributions from both the product of the two 0{L^) terms and the product of the 
0{L^) and 0{L'^) terms, we see the first deviation of Misner data from Brill-Lindquist data: Whereas these terms are 
non-zero for the Misner case, they vanish for the Brill-Lindquist case as the 0{L^) and 0{L^) terms are absent from 
the Weyl scalars for Brill-Lindquist initial data. Because the energy flux in gravitational waves is at 0{L^) — arising 
from the product of the 0{L^) terms — the energy flux in the Misner case deviates from its Brill-Lindquist counterpart 
at a relative order of 0{L^). Specifically, the relation of the energy fluxes in gravitational waves is given by 



^Misner J^^^ 



37r8 C2(3) \R 



(78) 



Note, that this relation holds separately for the flux in 'I'o and the flux in ^'4, and also for the total flux. 

Notice, that we can compare directly the radiation content of Misner data [Eqs. and ([7S)) ] and Brill-Lindquist 
data [Eqs. lpD|) and ([77]) ] using the BB scalar. We first notice that to 0{L^) (or, equivalently, to 0{zq)), the BB 
scalars are the same for both initial data sets. That is, the quadrupole radiations are the same. The reason for 
that is that our condition (|75p for the equality of the Coulomb curvatures is identical to the condition that the 
quadrupole amplitudes of the perturbations are equal, (zq/M) = 2-y/K2(Mo) [13 ■ Namely, by setting the Coulomb 
curvatures equal, we also set the quadrupole contributions to the BB scalar equal. We only see a difference at the 
next order: ^bl has a contribution at 0(2:0) (or, equivalently, at 0(L^)), whereas has a contribution at 0{z^^^) 
(or, equivalently, at 0{LP')). This implies that at small black-hole separations, or L <C Af, Brill-Lindquist initial data 
contain less gravitational radiation than Misner initial data. The same conclusion was also obtained by Abrahams 
and Price 13], who found that the Brill-Lindquist solution has a relatively smaller contribution because of higher 
multipole moments, as its geometry is more quadrupole dominated. Notably, because the difference between the two 
initial data sets enters at a lower order for the BB scalar than for the energy flux in gravitational waves, it would be 
easier to observe this difference numerically for the former than for the latter. The BB scalar is therefore a sensitive 
tool. Note, notwithstanding, that finding the BB scalar is not required for finding the Weyl scalars ^'4 and which 
is the objective. As we show in this paper, for the spacetimes we study here (which are both magnetic free) we can 
find the Weyl scalars. 
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